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Abstract 

We present a theoretical study of electrorotation (ER) of two spherical par- 
ticles under the action of a rotating electric field. When the two particles 
approach and finally touch, the mutual polarization interaction between the 
particles leads to a change in the dipole moment of the individual particle 

o. 

and hence the ER spectrum, as compared to that of the well-separated parti- 
cles. The mutual polarization effects are captured by the method of multiple 
images. From the theoretical analysis, we find that the mutual polarization 
j — , effects can change the characteristic frequency at which the maximum angular 

o 

velocity of electrorotation occurs. The numerical results can be understood 
in the spectral representation theory. 
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I. INTRODUCTION 

When a suspension of colloidal particles or biological cells is exposed to an external elec- 
tric field, the analysis of the frequency- dependent response yields valuable information on 
various processes, like the structural (Maxwell- Wagner) polarization effects |I]J2]. For most 
substances, the permittivity and conductivity are only constant over a limited frequency 
range. The general tendency is for permittivity to decrease, and conductivity to concomi- 
tantly increase, in a series of steps as frequency increases. These step changes are called 
dispersions. 

While the frequency- dependent response of biological cells can be investigated by the 
method of dielectric spectroscopy |§ , conventional dielectrophoresis and electrorotation (ER) 
analyze the frequency dependence of translations and rotations of single particles in an 
inhomogeneous and rotating external field, respectively [fU|. With the recent advent of 
experimental techniques such as automated video analysis || as well as light scattering 
methods |J, the motion of particles can be accurately monitored. In ER, the rotating field 
induces a cell dipole moment which rotates at the angular frequency of the external field. 
Any dispersion process causes a spatial phase shift between the induced dipole moment and 
the external field vector, giving rising to a torque which causes the rotation of individual 
cells [0. 

In the dilute limit in which a small volume fraction of particles or cells are suspended 
in a medium while the particles are well separated, we can concentrate on the ER response 
of individual particles by ignoring the mutual interactions between the particles. However, 
if the suspension is non-dilute, one may not ignore the interactions. Moreover, when the 
strength of the rotating electric field increases, the Brownian motion can be ignored and the 
polarized particles tend to aggregate along the rotating field even in the dilute limit. 

As an initial model, we will consider a pair of interacting particles dispersed in a sus- 
pension. When the two particles approach and finally touch, the mutual polarization inter- 
action between the particles leads to a change in the dipole moment of individual particles 



and hence the ER spectrum, as compared to that of the well-separated particles. We will 
present a theoretical study on the ER spectrum of two spherical particles in the presence 
of a rotating electric field, and capture the mutual polarization effects by the method of 
multiple images ||, recently developed to calculate the interparticle force between touching 
spherical particles. 

II. MULTIPLE IMAGE METHOD FOR A PAIR OF DIELECTRIC SPHERES 

We first consider an isolated spherical cell of complex permittivity e± = ei + ai/(iu) 
dispersed in a suspension medium of e 2 = e 2 + cr 2 /(io;), where to is the frequency of the 
external field Eq, and i = \J— 1. If one applies a rotating electric field, there is in general 
a phase difference between the induced dipole moment and the applied field vector, and 
thus a net torque is acting on the cell. The spin friction (analogous to the Stokes drag 
for translational motions) causes a steady-state rotational motion. In this case, the dipole 
moment of the isolated cell p is 

p = h 2 D 3 bE , (1) 

where b — (ei — e 2 )/(ii + 2e 2 ) is the dipole factor, and D the diameter of the cell. 

To account for the effect of multipole interaction on the elect rorot at ion, we consider a 
pair of touching spherical cells at a separation R suspended in a medium. We will calculate 
the effect of multipole interaction on the dipole moment. For a rotating field, the total 
dipole moment of one cell of the pair is given by 

p* = P T (cos 2 9) + P L (sm 2 9) = \{Pt + P L ), (2) 

where 9 is the angle between the dipole moment and the line joining the centers of the 
spheres. The average over all possible orientations of the cells gives a factor 1/2. In Eq.(2), 
Pt (Pl) is the transverse (longitudinal) dipole moment, being perpendicular (parallel) to 
the line joining the centers of the spheres. The transverse and longitudinal dipole moments 
are given by M 



where a satisfies the relation cosh a = R/D. For two touching particles, the dipole moment 
of a particle is given by 

f = h 2 D 3 b*E , (3) 

where the dipole factor of a pair is given by: 
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It is known that the angular velocity of electrorotation of a particle is related to the 
dipole factor (b or b*) as follows 

fi = —F(e2,r),E )Im[Dipole Factor] (5) 

where F is a function of e 2 , the dynamic viscosity 77 of the medium, as well as the electric- 
field magnitude E , and 7m[- • •] denotes the imaginary part of [•••]• For a spherical particle, 
F(e2,rj,E ) = e2El/2rj. Regarding r/, the spin friction expression suffices for an isolated 
particle. However, for two interacting particles, we must consider the more complicated 
suspension hydrodynamics. 

III. SPECTRAL REPRESENTATION AND DISPERSION SPECTRUM 

In a recent paper |jj , we studied the dielectric behavior of cell suspensions by employing 
the Bergman-Milton spectral representation of the effective dielectric constant |1(J. By 



means of the spectral representation, we derived the dielectric dispersion spectrum in terms 
of the electrical and structure parameters of the cell models. The spectral representation is a 
rigorous mathematical formalism of the effective dielectric constant of a composite material. 
The essence of the spectral representation is to define the following transformations. If we 
denote a complex material parameter 



'-if' (6) 

then the dipole factor b* admits the general form 

&* = E-^ (7) 

n ° °n 

where n is a positive integer, i.e., n — 1,2, ..., and F n and s n are the n-th microstructure 
parameters of the composite material JTH . 



Thus the spectral representation offers the advantage of the separation of material param- 
eters (namely the dielectric constant and conductivity) from the cell structure information, 
thus simplifying the study. From the spectral representation, one can readily derive the 
dielectric dispersion spectrum, with the dispersion strength as well as the characteristic fre- 
quency being explicitly expressed in terms of the structure parameters and the materials 
parameters of the cell suspension |J. The actual shape of the real and imaginary parts of 
the permittivity over the relaxation region can be uniquely determined by the Debye relax- 
ation spectrum, parametrized by the characteristic frequencies and the dispersion strengths. 
So, we can study the impact of these parameters on the dispersion spectrum directly. In 
what follows, we further express the dipole factor b and b* in the spectral representation. 
The dipole factor b admits the form b = Fi/(s — si) in the spectral representation, where 
Fi = —1/3 and Si = 1/3. For the dielectric dispersion spectrum, we further define two 
dimensionless real parameters || 

s = (1 - ei/e 2 ) _1 , and t = (1 - ai/o- 2 ) _1 , 

where s and t are the dielectric and conductivity contrasts respectively. After simple ma- 
nipulations, Eq.(7) becomes 

' Fl + ^%T (8) 



s - si 1 + if jf, 

where f c is the characteristic frequency, at which the maximum angular velocity of electro 
rotation occurs, and 5e± is the dispersion magnitude B: 
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Similarly, the dipole factor of a pair of interacting particles b* can be rewritten in the 
spectral representation as 

oo / rp(T) p(L) \ 
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where 
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In the above derivation, we have used the identity 
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We should remark that the spectral representation of b* [Eq.(lO)] is an exact transformation 

|n| of the multiple image expression [Eq.(4)]. Thus, the spectral representation of b* consists 



of a discrete set of simple poles (see Fig. 5 below). The m = 1 pole s\ (s\ ) for the 
longitudinal (transverse) field case deviates significantly from s = 1/3 while s[ becomes 
small in the touching limit, a — > 0. In what follows, we will show that this pole gives a 
significant contribution to the ER spectrum at low frequency. As m increases, however, the 
series s£P and sffl accummulate to s — 1/3, giving rise to a dominant contribution near the 
isolated-sphere frequency f c . Moreover, each term in the spectral representation expression 
of b* can be rewritten as, 

pCO jr(T) s at) 



(T) "" (T) ' -f/AT) 

S Sm S Sm -L ~r I J / Jmc 

p( L ) pW Se^ 



S Sm S Sm -L t %j j Jmc 

where 5effl (or Seffl) and f££) (or f^)) are the dispersion magnitudes and the characteristic 
frequencies of the transverse (longitudinal) field cases, obtained respectively by replacing 



Fi and si in the expressions of 5e x and f c in Eq.(9) with Fffl (or Fffl) and s^} (or s^). 
Note that / lc is significantly lower than f c . The ER spectrum of two spheres will consist 
of a series of sub-dispersions. In what follows, we will investigate the effect of multipole 
interaction on the electrorotation of cell, as shown in the next section. 

IV. NUMERICAL RESULTS 

For convenience, we let F(e2,t],E ) = 1 in our numerical calculation, both for the 
isolated-particle and touching-particle cases, in order to study the effect of multipole in- 
teraction on the ER spectrum. It is because we will focus on the frequency dependence 
of the ER spectrum. On the other hand, we do not know the dynamic viscosity for two 
approaching spheres. 

Fig.l is plotted to compare the two cases for three different conductivity contrast t. It 
is evident that, for both cases, a small \t\ yields a high characteristic frequency, at which 
the peak locates. Moreover, a smaller \t\ gives a larger peak value. Given a constant t, the 
touching-cell case predicts a shift of the peak location to lower frequency than that predicted 
by the isolated-cell case. At the same time, the peak for the touching-cell case broaden. This 
behavior arises from the effect of the multipole interaction. In addition, if t is given, the two 
cases also predict different behavior within the low-frequency region, but almost the same 
behavior within the high-frequency region. Especially for small \t\ (e.g., t = —1/900), a 
second peak occurs at lower frequency. This is in accord with the prediction of the spectral 
representation. 

Fig. 2 is plotted to exhibit the effect of the material parameter s on the two cases. It is 
observed that, for both cases, a larger s gives a larger peak value. Moreover, the peak of 
different s appears at almost the same frequency for both cases. In summary, s may affect 
the peak value, but not the peak location. 

Fig. 3 is plotted for three different values of the medium conductivity a 2 at s = 1.1 and 
t = —1/90. For a larger 02, we get a larger characteristic frequency at which the peak 
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occurs. It is evident that a-i may affect the peak location, but not the peak value. 

In Fig. 4, as R/D « 1 (e.g., R/D = 1.0333), the deviation between the two cases is 
evident. In other words, the multipole interaction does play an important role in the ER 
spectrum and the effect cannot be neglected for touching particles. However, as the separa- 
tion increases, say, R/D = 2, both cases predicts the same ER spectrum. From the results, 
we would say that the effect of multipole interaction may be neglected at R/D > 2. 

In Fig. 5, we plot the spectral representation for the two sets of poles s^ (or s%p) for 
m = 1 to 100. As m increases, sffl increases up to 1/3, while s£p decreases towards 1/3. For 
small R/D ratio (say, R/D = 1.0333), the longitudinal field plays an more important role 
in determining electrorotation spectrum than the transverse field. However, if R/D ratio is 
large enough, say, R/D > 2, (results not shown here), the effect of multipole interaction is 
negligible both for the longitudinal and transverse field cases and the results for an isolated 
cell recover. In summary, as the two cells approach and finally touch, one must take into 
account the effect of multipole interaction on the electrorotation of cells. 

DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. In this work, we attempted a theo- 
retical study of the electrorotation of two approaching spherical particles in the presence 
of a rotating electric field. When the two particles approach and finally touch, the mutual 
polarization interaction between the particles leads to a change in the dipole moment of 
individual particles and hence the ER spectrum, as compared to that of the well-separated 
particles. The mutual polarization effects are studied via the multiple image method. From 
the results, we find that the mutual polarization effects can change the characteristic fre- 
quency substantially. We should remarked that the multiple image method for two dielectric 
spheres is only approximate |12[| , but the approximation is quite good [§]. More accurate 
calculations based on bispherical coordinates can be attempted. However, we believe that 
similar conclusions can be obtained. 



8 



Moreover, when the volume fraction of the suspension is large, the rotating electric field 
leads to plate-like aggregation of particles in the plane of the applied field. In this case, the 
electrorotation spectrum can be modified further due to the local field effects arising from 
the many-particle system. To this end, a first-principles approach can be used to handle the 
many-particle and multipole interactions [pf]. 

Regarding the dynamic viscosity, the spin friction expression for an isolated spherical 
cell suffices. However, for many cells interacting in a suspension, we must consider the 
hydrodynamics. Regarding the dynamic viscosity for two touching spheres, we must consider 
the suspension hydrodynamics of two approaching spheres. This is a formidable task and 
should be a topic for future work. 
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FIGURES 
FIG. 1. The minus dipole factor (being proportional to the angular velocity of electrorotation) 

plotted versus frequency both for an isolated cell (dashed lines) and touching cells (solid lines) for 

three different conductivity contrasts t at s = 1.1 and 02 = 2.8 x 10 iSm -1 . 

FIG. 2. Similar to Fig.l but for three different dielectric contrasts s at t = —1/90 and 
02 = 2.8 x lO~ 4 Sm~ 1 for an isolated cell (lines) and touching cells (symbols). 

FIG. 3. Similar to Fig.l but for three different medium conductivities 02 at s = 1.1 and 
t = -1/90. 

FIG. 4. Similar to Fig.l but for three different separation ratios R/D at s = 1.1, t = —1/900 
and o- 2 = 2.8 x lO^Sm -1 . 

FIG. 5. The pole spectrum for two approaching spheres for several separation ratios R/D for 
the longitudinal (open symbols) and transverse (filled symbols) cases. Note that the longitudinal 
(transverse) part of the spectrum are on the small (large) s side of the spectrum. The lines are 
guides to the eyes. 
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Fig.2 
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Fig.3 
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